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Captures of stellar-mass compact objects (COs) by massive (~ 10^ M©) black holes (MBHs) 
are potentially an important source for LISA, the proposed space-based gravitational-wave (GW) 
detector. The orbits of the inspiraling COs are highly complicated; they can remain rather eccentric 
up until the final plunge, and display extreme versions of relativistic perihelion precession and Lense- 
Thirring precession of the orbital plane. The amplitudes of the strongest GW signals are expected 
to be roughly an order of magnitude smaller than LISA's instrumental noise, but in principle (i.e., 
with sufficient computing power) the GW signals can be disentangled from the noise by matched 
filtering. The associated template waveforms are not yet in hand, but theorists will very likely 
be able to provide them before LISA launches. Here we introduce a family of approximate (post- 
Newtonian) capture waveforms, given in (nearly) analytic form, for use in advancing LISA studies 
until more accurate versions are available. Our model waveforms include most of the key qualitative 
features of true waveforms, and cover the full space of capture-event parameters (including orbital 
eccentricity and the MBH's spin). Here we use our approximate waveforms to (i) estimate the 
relative contributions of different harmonics (of the orbital frequency) to the total signal-to-noise 
ratio, and (ii) estimate the accuracy with which LISA will be able to extract the physical parameters 
of the capture event from the measured waveform. For a typical source (a IOA^q CO captured by a 
10^ Mq MBH at a signal-to-noise ratio of 30), we find that LISA can determine the MBH and CO 
masses to within a fractional error of ~ 10~ , measure S/M^ (where S and M are the MBH's mass 
and spin) to within ~ lO"**, and determine the location to the source on the sky to within ~ 10"'^ 
stradians. 
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I. INTRODUCTION 

Captures of stellar-mass compact objects (COs) by massive (~ IO^Mq) black holes (MBHs) in galactic nuclei 
represent an important potential source for LISA, the proposed space-based gravitational- wave (GW) detector [1]. 
The capture orbits, which can remain rather eccentric right up to the final plunge, display extreme versions of 
relativistic perihelion precession and Lense-Thirring precession (i.e., precession of the orbital plane due to the spin of 
the MBH) , as well as orbital decay. Ryan [2] has illustrated how the measured waveforms can effectively map out the 
spacetime geometry close to the MBH. 

Rate estimates indicate that the strongest detectable sources will he D ^ 1 Gpc from us [3] , implying the measured 
GW strain amplitude h{t) will be roughly an order of magnitude smaller than the strain amplitude n{t) due to 
detector noise. Nevertheless, because the capture waveform will be "visible" to LISA for 10^ cycles, in principle 
(i.e., with infinite computing power) it should be possible to dig these signals out of the noise using matched filtering. 
While until now, theorists have not been able to calculate capture waveforms with the accuracy required for matched 
filtering, great progress has been made on this problem in recent years (see the recent review by Poisson [4] and 
references therein), and it seems very likely that sufficiently-accurate template waveforms will be available before 
LISA'S planned launch in ~ 2011. 

Nevertheless, there are many theoretical issues related to captures that must be addressed even before sufficiently- 
accurate templates are in hand, as the final LISA mission specs are being driven substantially by the requirement 
that LISA be sensitive enough to see captures (see the LISA Science Requirements Document at the LISA Sources 
and Data Analysis team website [5]). Among the most pressing questions is how to efRciently search through the data 
for capture signals. Given the large dimensionality of the space of capture waveforms, and the length and complexity 
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of the individual template signals, it is clear that straightforward matched filtering — using a huge set of templates 
that cover parameter space like a net — would require vastly more computational power than is practical. Instead, 
we need to develop suboptimal alternatives to coherent match filtering and to estimate the sensitivity of LISA when 
using these methods. 

This is the first in a series of papers designed to address such data analysis problems, for LISA capture sources. 
Here we introduce a family of relatively simple inspiral waveforms, given in nearly analjrtic form [i.e., up to solutions 
of ordinary differential equations (ODEs)], that should roughly approximate the true waveforms and that include their 
key qualitative features. Because our approximate waveforms are given in (nearly) analytic form, we can generate 
vast numbers of templates extremely fast — a feature that we expect to be crucial in performing Monte Carlo studies 
of search techniques. 

Our capture waveforms are based on the lowest-order, quadrupolar waveforms for eccentric-orbit binaries derived 

by Peters and Matthews [6] , but the orbits are corrected to include the effects of pericenter precession, Lcnsc-Thirring 
precession, and inspiral from radiation reaction [all calculated using post-Newtonian (PN) formulae] . Our waveforms 
have the right dimensionality (except that in practice we will neglect the spin of the CO) and the same qualitative 
features as the true waveforms, except (i) they do not exhibit the extreme "zoom-whirl" behavior of true fully- 
relativistic waveforms just prior to plunge [7,8], and (ii) they approximate as constant the angle A between the CO's 
orbital angular momentum and the MBH's spin, whereas the correct evolution would show a small secular change in 
A. However, we do not expect these two (missing) effects to be very important for our purpose, since, again, the latter 
is small (at least for circular orbits [9]), while extreme zoom- whirl behavior occurs only for rather eccentric orbits 
extremely close to plunge. [The number of "whirls" per "zoom" grows only as — ln(p — Ps) as the orbit's semi-latus 
rectum p approaches the plunge value Ps [7,8], and so stretches of waveform exhibiting extreme zoom- whirl behavior 
will likely contribute only a very small fraction of the total signal-to- noise ratio (SNR).] 

In this paper, we use our approximate waveforms to take a first cut at estimating the accuracy with which LISA 
should be able to extract the capture source's physical parameters, including its distance and location on the sky, the 
masses of both bodies, and the spin of the MBH. We also illustrate in detail the relative contributions of different har- 
monics (of the orbital frequency) to the total SNR of the waveform. For these calculations, besides using approximate 
waveforms, we also use a low- frequency approximation to the LISA response function and an approximate version of 
the Fisher-matrix inner product. Since our methods are so approximate, the results should be considered illustrative 
rather than definitive. 

The plan of this paper is as follows. In Sec. II, to put the present work in context, we briefly summarize the 
previous literature on this topic and indicate the several areas of active theoretical research. In Sec. Ill wc present 
our approximate, (nearly-) analytic waveforms. Our scheme for incorporating the LISA response function and the 
noise closely follows Cutler [10], but we give enough detail that the reader could easily employ the response function 
detailed in Cornish and Rubbo [11] (which is more accurate at high frequencies). Sec. IV provides a brief review of 
signal analysis, partly to explain our conventions. In Sec. V we display plots showing the relative contributions of 
different harmonics to the total SNR. Wc show that, even for relatively modest final eccentricity, the higher harmonics 
contribute significantly. Finally, in Sec. VI, we present estimates of how accurately LISA can determine the physical 
parameters of capture systems. We emphasize that our treatment is highly modular, allowing for simple improvements 
of the various approximations. 

A few details of our analysis are left to the Appendices. In Appendix A we derive a simple expression for the 
contribution to the pericenter precession from the spin of the MBH, and show its equivalence to the standard result. 
Appendix B compares the magnitudes of the various PN terms in our orbital evolution equations. Since the near- 
plunge capture orbits are highly relativistic, higher-order terms are comparable in magnitude to the lower-order ones, 
as one would expect. Estimates of the magnitudes of effects related to the CO's spin are given in Appendix C. The 
spin of the CO may have marginally important effects on the templates (for rapidly rotating COs), but for simplicity 
we leave these out of the rest of our analysis. (They could be put back in rather easily.) 

In later papers will turn to the problem that is the main motivation for this work: designing a practical algorithm 
for digging capture waveforms out of the LISA noise. There we will use our approximate waveforms to estimate the 
scheme's sensitivity, compared to an optimal search with infinite computing power. 

Our index notation is the following. Indices for vectors and tensors on parameter space are chosen from the 
beginning of the Latin alphabet (a, b.c, . . .). Vectors and tensors on three-dimensional space have indices chosen from 
the middle of the Latin alphabet {i.j, k, . . .), and run over 1, 2, 3; their indices are raised and lowered with the flat 
3-metric, rjij. We use Greek indices {a,f3, . . .), running only over I, II, to label the two independent gravitational 
waveforms that LISA effectively generates. (No four-dimensional, spacetime indices occur in this paper.) 

Throughout this paper we use units in which G = c=l. 
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II. SUMMARY OF PREVIOUS AND ONGOING THEORETICAL WORK 



Most, if not all, nucleated galaxies harbor MBHs in their centers [12,13]. The MBH's gravity dominates the local 
stellar dynamics within a cusp radius Tc = M/cr^, where M is the MBH mass and Cc is the onc-dimcnsional velocity 
dispersion of stars inside the cusp. A typical MBH with M = IQ^Mq would have Tc ~ 1 pc. The total mass of stars 
inside the cusp is typically of order the MBH mass [3]. Captures occur when two stars in the cusp undergo a close 
encounter, sending one of them into the "loss cone." These are orbits that pass sufficiently close to the MBH that the 
timescale on which the CO tends to spiral into the MBH due to gravitational radiation reaction is shorter than the 
timescale on which the CO is scattered back out by other stars. 

Because EISA's sensitivity band is centered at / ~ 3 x 10~^Hz, the MBHs most "visible" to LISA are those with 
mass M ~ IQ^Mq. To avoid tidal disruption, while being close enough to the MBH to emit GWs in the LISA 
frequency band, the captured star must be either a white dwarf (WD), neutron star (NS), black hole (BH), or a very 
low mass main-sequence star (LMMS) [3]. 

Early estimates of capture rates and LISA SNRs were made by Hils and Bender [14], who considered the capture 
of IMq objects. More recent rate estimates [3] suggest that, while the total capture rate is dominated by LMMSs 
and WDs, LISA's detection rate should be dominated by captures of ~ IOMq BHs. This is partly because the BHs, 
being more massive, can be "seen" to greater distance, and partly because two-body relaxation enhances the density 
of BHs nearer the MBH [15]. (Two-body stellar collisions tend to equalize kinetic energies, causing heavier stars to 
sink to the center of the cusp.) 

The first extended look at data analysis for capture sources was taken by Finn and Thorne [16]. They simplified 
the problem by restricting to the case of circular, equatorial orbits, but for this case they were able to calculate 
the correct relativistic orbits and waveforms, and they showed how the LISA SNR accumulates over the last year 
of inspiral — during which typically ^ 10'"' GW cycles arc emitted for a range of CO masses and MBH masses and 
spins. Their plots illustrate the salient fact that, typically, the entire last year of inspiral contributes significantly to 
the SNR. This is because, one year before plunge, the CO is already quite close to the MBH. Indeed, we shall see 
below that considerable signal-to-noise can accumulate even before the final year. 

A more realistic treatment of the capture problem must incorporate the facts that (i) capture orbits will generally 
be non-equatorial (i.e., the CO's orbital angular momentum will not be aligned with the spin of the MBH), and (ii) 
a fair fraction of the inspiraling orbits will remain moderately eccentric right up until the final plunge. The latter 
fact may seem surprising, since it is well known that gravitational radiation reaction tends to circularize orbits rather 
efficiently. "'^ The point, however, is that when the COs enter the loss cone, their orbits are initially extremely eccentric: 
1 — Cinit ^ 10~^ — 10~^, typically, while the initial pericenter distance is only rp^init ^ 8 — lOOM [18]. Given the CO's 
initial trajectory, just after scattering into the loss cone, we would like to calculate the eccentricity at the last stable 
orbit, clso- For non-spinning MBHs, at least, this is straightforward. We find that e^so > 0.1 if rp, init < 20. OM, 
^LSO > 0.2 for Tp^init <i 12. 8M, and Clso > 0.3 for r^^init ^ 9.2M. [These estimates were obtained as follows. In 
the test particle limit, let ri and r2 be the turning points (pericenter and apocenter) of the radial motion, where r is 
the standard radial coordinate in Schwarzschild. Define p and e by ri(= rp) = p/{l + e) and r2 = p/{l — e). Plunge 
occurs at p/M = 6 -|- 2eLSO [7]. Then pinit is given by 



(Of course, the upper limit in the integral should actually be slightly less than 1.0 — say, e = 0.99995 — but since 
the integrand is smooth as e — > 1.0, it makes no practical difference if we simply approximate the upper limit as 
1.0.) The derivative dp/de = p/e due to radiation reaction was calculated numerically by Cutler, Kennefick, and 
Poisson [7] for orbits near the horizon. We used the results from Fig. 1 of [7] to integrate (roughly, using a pencil and 
ruler) dp/de backwards (in time) from plmige to p/M = 12. We then used the lowest-order post-Newtonian result [6] 
dp/de = (12/19)(]3/e)(l + |e^)/(l + ^^jc^) to continue the integration backwards to e = 1.0.] 

Based on Freitag's Monte Carlo simulation of capture events in our Galaxy (Fig. 1 of [18]), we then estimate 
that roughly half the captures of ~ IOMq BHs (which, again, should dominate LISA's detection rate) should have 
eLSO ^ 0.2. Note that a year or two before the final plunge the eccentricity of such captures will, in fact, be 
significantly larger than blso (as illustrated in Figs. 2,3 below). 



^Except very close to plunge, where the very strong-field potential tends to decrease the rate of circularization, and may even 
reverse the sign of de/dt — of. [7,17]. 




(1) 
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Capture sources are unlike some other LISA sources (e.g., galactic WD- WD binaries or MBH-MBH binaries at 
high redshift), in that they may lie quite near the margin of detectability, given LISA's current design specifications. 
(Put another way, a modest change in the height or location of the noise floor may determine whether or not these 
sources are detected.) Since detecting capture sources is very high priority for LISA, it has been a high priority for 
LISA'S Sources and Data Analysis team ["Working Group 1" (WGl)] to make as much progress as possible studying 
capture sources before finalizing the LISA design [19]. This motivates current research on several fronts, including 
work to (i) improve estimates for event rates and for the distribution of source parameters (especially masses and 
initial pericenter distances); (ii) solve the radiation reaction problem to determine the true orbit, and construct the 
corresponding waveforms; (iii) investigate what science can be done with these sources (both astrophysics and tests of 
fundamental physics); (iv) understand the limits on capture detection due to "source confusion", i.e., the background 
"noise" caused by other, unresolved capture sources; and (v) construct strategies to dig the capture waveforms out of 
the instrumental and confusion noise. The present work addresses issues (iii) and (iv) above, while later papers will 
address problem (v). 

Parameter estimation with LISA [clearly bearing on above issue (iii)] has been looked at systematically for WD- 
WD binaries by Peterseim et al. [20] and Cutler [10], and for mergers of MBH pairs by Cutler [10] and Vecchio [21]. 
No comparable analysis has yet been done for capture sources. For captures, some initial estimates of parameter 
estimation accuracy were made by Poisson [22] and Ryan [2] (the latter's main interest being to test alternative 
gravitation theories). However, both Poisson and Ryan used extremely simplifying approximations: they both took the 
inspiral orbits to be circular and equatorial a priori (effectively reducing the number of unknown system parameters, 
while leaving uninvestigated the significance of perihelion precession and Lense-Thirring precession for parameter 
extraction), and they did not incorporate in their signal models the amplitude and phase modulations that arise from 
LISA'S orbital motion (which LISA will use to determine the source position). By comparison, our treatment is far 
more realistic. While our results are also approximate, we believe they should at least give correct order-of-magnitude 
estimates of LISA's parameter estimation accuracy (while it seems doubtful that the earlier estimates can be trusted 
even at that level). 

III. OUR MODEL INSPIRAL WAVEFORM 

We approximate the CO-MBH system as being, at any instant, a Newtonian-orbit binary emitting a Peters- 
Matthews (i.e., lowest order, quadrupolar) waveform. We then use post-Newtonian (PN) equations to secularly 
evolve the parameters of the orbit. In particular, we include orbital decay from radiation reaction, pericenter preces- 
sion, and Lense-Thirring precession of the orbital plane. The modulation of the waveform's amplitude and phase due 
to Lense-Thirring precession has been described (in the context of circular-orbit binaries) by Apostolatos et al. [23] . 
The motion of the LISA detector introduces additional modulations; our handling of these closely follows that of 
Cutler [10]. Cutler's treatment does not account for the decrease in LISA's sensitivity at frequencies / ^ 10 mHz 
(where the GW wavelength is smaller than the detector's armlength) an effect recently accounted for by Cornish and 
Rubbo [11]. It would be straightforward to repeat our analysis using the Cornish and Rubbo formalism. However, 
since most of the SNR for astophysically-relevant capture sources will accumulate at frequencies below 10 mHz, we 
expect our high-frequency approximation will have only a modest impact on the results. 

We emphasize that our treatment is highly modular. E.g., while our choice of physical variables is a slight gen- 
eralization of the variables used in [23] and [10], it would be straightforward to re- write our waveforms using a 
parameterization along the lines of Buonanno, Chen, and Vallisneri [24] (hereafter, BCV), who found a particularly 
convenient way of parametrizing circular-orbit binaries with spin. (It seems the BCV parametrization should be read- 
ily extendible to eccentric orbits.) Similarly, it would be straightforward for us to treat the LISA noise and response 
function along the lines of Cornish and Rubbo [11]. We do not implement either of these treatments here because they 
did not seem of critical importance, and because our work was already well underway when they became available. 

A. Principal axes 

In this paper, we adopt a mixed notation for spatial vectors, sometimes labelling them with spatial indices 
fc, . . .), but sometimes suppressing the indices and instead using the standard 3 — d vector notation: an over-arrow 
(as in A) to represent a vector. A- B to represent a scalar ("dot") product, and Ax B to represent the vector ("cross") 
product. An over-hat (as in h) will indicate that a vector is normalized, i.e., has unit length. We trust our meaning 
will always be clear, despite this mixed notation. 
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Let h be the unit vector pointing from the detector to the source, and let L{t) be the unit vector along the CO's 
orbital angular momentum. We find it convenient to work in a (time-varying) wave frame defined with respect to fi 
and L{t). We define unit vectors p and q by 

j5 = (n X L)/\h X Z/|, 

q = pxfi, (2) 
based on which we then define the two polarization basis tensors 

H^^j{t)=Mj+mPr (3) 
The general GW strain field at the detector can then be written as 

h,j{t)=A+{t)H+{t) + A-{t)H^.{t), (4) 
where A'^{t) and {t) are the amplitudes of the two polarizations. 

B. Peters-Mathews waveforms 

In the quadrupole approximation, the metric perturbation far from the source is given (in the "transverse/traceless" 
gauge) by [25] 

hij = {2/D) {PikPji - IPijPki)!"' (5) 

where D is the distance to the source, the projection operator Pij is given by Pij = rjij — hifij , and is the second 

time derivative of the inertia tensor. In this paper we work in the limit of small mass ratio, fi/M <C 1, where fi and M 
are the masses of the CO and MBH, respectively. In this limit, the inertia tensor is just P^{t) = /i {t)r^ {t) , where r 
is the position vector of the CO with respect to the MBH. 

Consider now a CO-MBH system described as a Newtonian binary, with semi-major axis a, eccentricity e, and 
orbital frequency v = (27rM)~^(M/a)"^/^. Let ei and 62 be orthonormal vectors pointing along the major and minor 
axes of the orbital ellipse, respectively. Since the orbit is planar, P^ has only 3 independent components: P^, P^ , and 
P"^, and as the motion is periodic, we can express P^ as a sum of harmonics of the orbital frequency v. P^ = . 

We next denote 

b =/i2 

Cn^l{l'n'+I^'). (6) 

Peters and Matthews showed [6] that 

a-n = -nA[Jn-2{ne) - 2eJ„_i(ne) -I- (2/n) J„(ne) + 2eJ„+i(ne) - J„+2(ne)] cos[n'I>(t)], 
hn = -nA{l - e^fl'^ [J„_2(rie) - 2 J„(ne) J„+2(ne)] sin[n$(t)], 

c„ = 1AJn{ne) cos[n$(t)], (7) 



where 



A = (27rzyM)2/3^, (8) 



J„ are Bessel functions of the first kind, and $(t) is the mean anomaly (measured from pericenter). For a strictly 
Newtonian binary we have 

$(t) = 2T:v(t - to) + $0, (9) 



5 



where $o is the mean anomaly at to- Decomposing Eq. (4) into n-harmonic contributions and using Eq. (5), one then 
easily obtains explicit expressions for the n-harmonic components of the two polarization coefficients, = 
and A~ = A~ . They are 

A+ = -[l + iL- [a„ cos(27) - 6„ sin(27)] + [I - [L ■ fifjcn, 

= 2{L ■ n) [bn cos(27) + a„ sin(27)] , (10) 

where 7 is an azimuthal angle measuring the direction of pericenter with respect to a; = [— n + I/(L-n)]/[l — (L-n)^]^/^. 



C. LISA'S response function 

With its three arms. LISA hmctions as a pair of two-arm detectors, outputting two orthogonal signals. Let l\,l2: l\ 
be unit vectors, each along one of EISA's three arms, and let L be EISA's average arm length. Let also Li{t) be the 
length of the i'th arm when LISA measures an incident GW, and denote 5Li{t) = Li{t) — L. We refer to the two-arm 
detector formed by arms 1 and 2 as 'detector E' The strain amplitude in this detector is given by 

hi{t) = [5Li{t) - SL2it)] /L = ^ hij{t){lili - (11) 

The second, orthogonal signal is then given by [10] 

hn{t) = 3-^/' [SL,{t) + SL2{t) - 2SLs{t)] = ^ hij{t){l{l{ + l^i - 2lili). (12) 

For GW wavelengths much larger than the LISA arm length, hi{t) and hii{t) coincide with the two "Michelson 
variables" [26], describing the responses of a pair of two-arm/90° detectors. We can then write hi{t) and hii{t) as a 
sum over n-harmonic contributions, 

/i„(t) = {a = 1,11), (13) 

n 

where 

hc^At) = [F+it)A+{t)+F-{t)A:it)] . (14) 

Here, A^'^ {t) arc the two polarization coefficients [given, in our model, by Eq. (10) above], the factor \/3/2 accounts 
for the fact that the actual angle between LISA arms is 60° rather than 90°, and F^'^ are the "antenna pattern" 
functions, reading [23,27] 

F+ = ^(1 + cos^ d) cos(2(/)) cos(2V') - cos6'sin(20) sin(2V'), 

F/ = ^ (1 + cos^ 6) cos{2(j)) sm{2tp) + cos 9 sin(2(;i) 008(2^)), (15a) 



F+ = 2 (1 + cos^ ^) sin(2^) cos(2V') + cos6'cos(20) sin(2V'), 

= ^(1 + cos^ 0) sin(2^) sin(2V') - cos6»cos(2(/)) cos(2^). (15b) 

In those expressions, (0, (p) is the source's sky location in a detector-based coordinate system and ?/' is the "polarization 
angle" describing the orientation of the "apparent ellipse" drawn by the projection of the orbit on the sky — see Fig. 
1 in Ref. [23] and the explicit relation (17) given below. 

It is more convenient to express the above response function in terms of angles defined not in the rotating, detector- 
based system, but rather in a fixed, ecliptic-based coordinate system. The angles 9, are related to 9s, (ps — the source 
location in an ecliptic-based system — through 
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1 \/3 

cos6'(t) = - cos 61s - ^ sines cos[<^o + 27r(t/T) - 0s], 



(j){t) =ao + 2'!r{t/T) + tan"^ 



V3 cos 6*5 + sin 6*5 cos[(^o + 27r(t/T) - (/)s] 
2 sin 9s sin[^o + 27r(Vr) - (f>s] 



(16) 



where T = 1 year and (^o^ao are constant angles specifying, respectively, the orbital and rotational phase of the 
detector at t = 0. (See Cutler [10] for a complete definition of these angles; note, though, that the angle ao = in 
this paper is referred to as ao = in Cutler [10].) 

Next, we express the polarization angle tp in terms of 6s, (f)s and 6l, (t>L — the direction of the CO's orbital angular 
momentum, L{t), in the ecliptic-based system. We have 



I 1 %/3 

tanV' ~ ^ 2 ^'^^^i^ ^ sin^L cos[(^o + ^^{t/T) — — cosO{t) [cos^l cos0s + sin^x sin^s cos((/)i — (j)s) 

1 \/3 

- sihOl sin 6*5 sin(0i — (ps) ^ cos((/)o + 27rf/T) {cos0l sin 6*5 sin (ps — cos 6*5 sihOl sin^^) 

- sin(0o + 27ri / T) (cos 63 sin 6l cos (pL — cos 6 l sin 63 cos (ps ) 



I 



(17) 



For concrctcncss wo shall hereafter take (pQ = oiq = 0, but one could specify any other value as appropriate. 

Note that the angles 9l, <t>L are not constant, since L precesses about the MBH's spin direction S. Let 9k, 4'k be 
the direction of S in the ecliptic-based system ('K' standing for 'Kerr'); let also A be the angle between L and S, and 
a{t) be an azimuthal angle (in the orbital plane) that measures the precession of L around S: Specifically, let 



L = S cos A 



Scoti9K . . 

■ sm A cos a 



Sxz 



sin A sin a, 



(18) 



sin 9 k sin 9 k 

where i is a unit vector normal to the ecliptic. Then the angles ^i(f), (f>L{t) are given in terms of 9k, ^k, A, a{t) as 

cos (t ) = cos 9k cos a + sin 9k sin A cos a{t) , 

sin 9l (t) cos (pi (t) = sin 9k cos (pK cos A — cos (pK cos 9k sin A cos a{t) + sin (pK sin A sin a{t) , 

sin 9 L{t) sin (pL{t) = sin 9k sin (pK cos X — sin (pK cos 9k sin X cos a{t) — cos(/>if sinAsina(t). (19) 



D. The pericenter angle 7 

As mentioned above, the angle 7 that appears in Eqs. (10) measures the direction of pericenter with respect to 
x=[-h + L{L ■ n)]/[l - {L ■ n)^] V^- With this definition, 7 is neither purely extrinsic nor purely intrinsic. (In the 
terminology of BCV, "intrinsic" parameters describe the system without reference to the location or orientation of 
the observer.) We will find it convenient to introduce a somewhat diflFerent convention for the zero-point of this angle: 
We shall define 7 to be the direction of pericenter with respect to L x S. Then 7 is a purely intrinsic quantity. 

Clearly, 7 and 7 are related by 

1 = 7 + 13, (20) 
where (3 is the angle from x oc [L{L • n) — n] to {L x S). It is straightforward to show that (3 is given by 

cos XL ■ h — S ■ h 



sin (3 ■ 



sinA[l- (i-n)2] 



1/2 ' 



„ n ■ is X L) , 

cos/J= , . ,1/2 - (21) 

sinA[l- (i-n)2] 

To evaluate /3(f) in practice, we shall need the following relations: 

S ■ h = cos 9s cos 6k + sin 6s sin 9k cos{(ps — (pK ) , (22) 
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, /A fx .^./, ,x.> S ■ ticosOk — cosOs . , . 

n- [b X L) = sm 6s sm[(pK — <ps) sin A cos a H : — sin A sin a, (23) 

sm ux 

and 

f - a - X , cos6's - S* • ncos6'K (S x z) ■ h . 

L ■ n = b ■ n cos A H : — sin A cos a H : — sin A sm a, (24) 

sm Ok sm 6k 

or, equivalently, 

L- n = cos6s cos6l + sm6s sin^j, cos{(f)s — (t>L)- (25) 

Note that the time-variation of • n is very smah in the extreme mass-ratio ease considered here: this quantity is 
constant to better than ~ {fi/M){S/M'^) (see Appendix C). In our model we shall approximate S — and hence S-ri — as 
strictly constant. 



E. Parameter space 

The two-body system is described by 17 parameters. The spin of the CO can be marginally relevant (see Appendix 
C), but in this paper we shall ignore it, leaving us with 14 parameters. We shall denote a vector in the 14-d parameter 
space by A" (a = 0, . . . , 13). We choose our parameters as follows: 

A« = (A°, . . . , A^^) = 

[to Inyu, InM, cq, 7o, $o, Ms = cost's, 4>s-, cos A, ao,iJ.K = cos6k, ^k, InW^)] • (26) 

Here, to is a time parameter that allows us to specify "when" the inspiral occurs we shall generally choose to to be the 
instant of time when the (radial) orbital frequency sweeps through some fiducial value fo (typically, we shall choose 
vq of order 1 mHz), /i and M are the masses of the CO and MBH, respectively, and S is the magnitude of the MBH's 
spin angular momentmn (so < S/M"^ !i !)• The parameters eo, 70, and $0 describe, respectively, the eccentricity, 
the direction of the pericenter within the orbital plane, and the mean anomaly — all at time to- More specifically, we 
take 7o to be the angle (in the plane of the orbit) from Lx S to pericenter, and, as usual, $0 to be the mean anomaly 
with respect to pericenter passage. The parameter ao = oi(t = to) [where a{t) is defined in Eq. (18)] describes the 
direction of L around S at to- The angles {63, (ps) are the direction to the source, in ecliptic-based coordinates; 
{(^K, ^k) represent the direction S of the MBH's spin (approximated as constant) in ecliptic-based coordinates; and 
A is the angle between L and S (also approximated as constant^). Finally, D is the distance to the source. 

The various parameters and their meaning are summarized in Table I. Fig. 1 illustrates the various angles involved 
in our parameterization. 

Note for simplicity we arc treating the background spacetime as Minkowski space, not Robertson-Walker. To correct 
this, for a source at redshift z, requires only the simple translation: M M(l + z), /j. ^ /x(l + z), S ^ S{1 + z)^, 
D ^ Dl, where Dl is the "luminosity distance" [28]. 

The parameters can be divided into "intrinsic" and "extrinsic" parameters, following BCV. Extrinsic parameters 
refer to the observer's position or orientation, or to the zero-of-time on the observer's watch. There are seven extrinsic 
parameters: the four parameters to, Us, <ps, and D correspond to the freedom to translate the same binary in space 
and time, and the three parameters fXK, (pK, and ao are basically Euler angles that specify the orientation of the 
orbit with respect to the observer (at to)- The intrinsic parameters are the ones that control the detailed dynamical 
evolution of the system, without reference to the observer's location or orientation. In our parametrization, the seven 
intrinsic parameters are In/x, InM, S/M"^, cos A, eo, 70, and $o- BCV observed (in the context of circular-orbit 
binaries with spin) that extrinsic parameters are generally much "cheaper" to search over than intrinsic parameters. 
We shall make good use of this important observation in further papers. 



^In reality, radiation reaction will impose a small time variation in A; however, this variation is known to be very small (See 
Ref. [9]) and we shall ignore it here. When a model of the time- variation of A is eventually at hand, it would be trivial to 
generalize our treatment to incorporate it: one would just need an equation for dX/dt, and in the parameter list A would be 
replaced by Ao — the value of A at time to- 
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FIG. 1. The MBH-CO system: setup and notation. AI and /i are the masses of the MBH and the CO, respectively. The axes labeled 
x—y—z represent a Cartesian system based on ecliptic coordinates (the Earth's motion around the Sun is in the x— y plane). The spin 
S of the MBH is parametrized by its magnitude S and the two angular coordinates d]f^<t>K^ defined (in the standard manner) based on 
the system x—y—z. Lit) represents the (time-varying) orbital angular momentum; its direction is parametrized by the (constant) angle A 
between L and 5, and by an azimuthal angle a{t) (not shown in the figure). The angle 7(4) is the (intrinsic) direction of pericenter, as 
measured with respect to L X S. Finally, 3>(t) denotes the mean anomaly of the orbit, i.e., the average orbital phase with respect to the 
direction of pericenter. 



F. Orbital evolution equations 

We evolve v{t), l{t), e(t), and a{t) using the following PN formulae: 

- = 2.., (27) 

^ = ii(^/M3)(2^Af:.)"/3(l - e2)-9/2{ [1 -t- (73/24)e2 + {37/96)e^] (1 - e^) 

+ {2-KMvf'^ [(1273/336) - (2561/224)e2 - (3885/128)6* - (13147/5376)6^] 
-{2^Mv){S/M'') cos A(l - e2)-i/2 [(73/12) + (1211/24)6^ 

+ (3143/96)6"*+ (65/64)6^]}, (28) 

-12'kvcos\{S/M'^){2ttMv){1 - e^)"^/^ (29) 
" (a*/M^)(1 - e^)-'^'^{2TTMvf/^ [(304 + 12l62)(l - e^){l + l2{2TTMvf''^) 

--^(27rMzy)2/3((8)(16705) + (12)(9082)62 ~ 252116-*)] 

+e{^l/M'^){S/M'^) cos A (27rAfiy)"/3(i „ ^^yi [(1364/5) + (5032/15)6^ + (263/10)e*] , (30) 
= Ativ{S/M^){2t:Mv){1 - e^)-^'"^. (31) 



de 6 
'di ~ ^15 



1 + i(27r;/M)2/3(l - e'^)-\26 - 156^) 



da 



Equations (28), (29), and (30) are from Junker and Schafer [29], except (i) the second line of Eq. (29) is from 
Brumberg [30] (cf. our Appendix A), and the last term in Eq. (28) — the term cx S/JVP — is from Ryan [31]. Eq. (31) is 
from Barker and O'Connell [32]. The dissipative terms dv/dt and de/dt are given accurately through 3.5PN order (i.e., 
one order higher than 2.5PN order, where radiation reaction first becomes manifest)."^ The non-dissipative equations, 



^The currently undetermined term in the 3.5PN expressions [see [33], in particular Eqs. (12)-(14) therein] does not show up 
in our calculation, since here we are ignoring terms that are higher-order in the mass ratio /i/M. 
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Ai 


In ^ 


Hn ofl ClO's TTlrlSR 


A2 


InM 


Qn ofl MBH's mass 




s/m'^ 


mil 0"n ifnrlp of I c:npr'i'fip 1 "^nin Jino'Tilar mmnPTli"mn of A/TRTT 


A4 


"^U 


6(^0), wlierG c(^) is the orbital 6CC6ntricity 


A^ 


/U 


'/(^o)^ where ^(t) is the angle (in orbital plane) between L x S and pericenter 


A6 


$0 


^(to), where ^(t) is the mean anomaly 


A^ 


^5 = COS ^5 


(cosine of) the source direction's polar angle 


A8 




azimuthal direction to source 


A^ 


COS A 


L ■ S{= const) 


Aio 


ao 


a{to), where a{t) is the azimuthal direction of L (in the orbital plane) 


A" 


= cos 6 k 


(cosine of) the polar angle of MBH's spin 


Ai2 




azimuthal direction of MBH's spin 


Ai3 




(In of) CO's mass divided by distance to source 



TABLE I. Summary of physical parameters and their meaning. The angles (Osjips) and {6K,<t>K) are associated with a spherical 
coordinate system attached to the ecliptic. L and 5 axe unit vectors in the directions of the orbital angular momentum and the MBH's 
spin, respectively. For further details see figure 1 and the description in the text. 

for drijdt and da/dt, are accurate through 2PN order.^ 

In solving the above time-evolution equations, the initial values (at time to) of v, 7, e, and a are just the 
parameters $0) I'O) 70) eo, and ao- 

We emphasize again that our treatment is highly modular: The PN expressions in Eqs. (28)-(31) could be replaced 
with improved ones as soon as higher-order PN expressions are available. Also, one might wish to improve these 
evolution equations using values from look-up tables, or results from numerical studies of the orbital evolution in Kerr 
(such as in [9]). 

G. Doppler phase modulation 

Doppler phase modulation due to LISA's orbital motion becomes important for integration times longer than a few 
weeks. We incorporate this effect by shifting the phase $(t), according to 

$(t) ^ $(t) + $^(t), (32) 

where 

i>^(t) = 27ri/(t)i?sin6iscos[27r(Vr) - 0s]. (33) 

Here i? = 1 AU = 499.00478 sec. 

H. Putting the pieces together 

The algorithm for constructing our approximate waveform is then: Fix some fiducial frequency and choose 
waveform parameters (to, In^, InM, S/M"^, eo, 70, $0, coaOs^ (j)Si cos A, uq^cosOk, 4'k,D). Solve the ODEs (27)- 
(31) for $(t), v{t), ^(t), e{t), a{t). Calculate 9L{t),4'L{t) using Eqs. (19) and then obtain ■)/'(<) from (17). Calculate 
j{t) from 7(t) using Eqs. (20) and (21). Use e{t) and i^{t) to calculate a„(t), 6„(t), c„(i) in Eqs. (7), remembering 
to include the Doppler modulation via n^{t) n[^{t) + $^(t)], a la Eqs. (32) and (33). Calculate the amplitude 



*In fact, the equations for d^/dt and da/dt are missing terms proportional to (iS/M^)^, which, according to usual "order 
counting" are classified as 2PN. However, this usual counting is misleading when the central object is a spinning BH: Because 
BHs are ultracompact, their spins are smaller than suggested by the usual counting, and the missing terms oc (5/M^)^ have, 
in fact, the same magnitude as 3PN terms. Similarly, the terms oc (iS/M^) in Eqs. (29) and (31) can be viewed as effectively 
1.5PN terms. 
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coefficients and the antenna pattern functions F^'^ using Eqs. (10) and (15), respectively. Then finally calculate 
ha{t) (for a = I, II) using Eqs. (14) and (13). 

Note that, in our treatment, pericenter precession and Lense-Thirring precession have no effect on the intrinsic 
signal (the signal in a frame that rotates with the system), since we always use the Peters-Matthews "lowest-order" 
waveforms. The effect of these motions is simply to rotate the binary system with respect to the detector. This relative 
rotation modulates the polarization angle ^ [which appears in the response functions F^'^' (t)] since depends on L, 
and it affects the amplitudes since the latter depend on L ■ n and 7. 

Figures 13-17 show some sample orbits and waveforms, obtained from employing the above algorithm. Figures 2 
and 3 demonstrate the evolution of orbits in the eccentricity-frequency plane. 




eccentricity 

FIG. 2. Evolution of orbits in our model, for a system composed of a IMq CO inspiralling into a IO^Mq (non-spinning) 
MBH. The dashed line represents the last stable orbit (LSO). Each of the solid lines shows the u — e trajectory of a system 
with given initial data (the orbit evolves in time "from bottom to top"). The four dots plotted along each trajectory indicate, 
from bottom to top, the state of the system 10, 5, 2, and 1 years before the LSO. 




0' ' ' ' ' ' ' ' ' ' 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

eccentricity 

FIG. 3. Same as in Fig. 2, for a WMq CO inspiralling into a IO^'Mq MBH. 
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IV. FORMALISM OF SIGNAL ANALYSIS 



This section briefly reviews the basic formulae of signal analysis, with application to LISA. We follow closely the 

treatment of Cutler [10]. In particular, our analysis is strictly valid only in the low- frequency regime, where the 
light-travel-time up and down one arm is much less than the gravitational wave period. (However, we expect it to be 
a reasonable approximation at higher frequencies as well). 

As discussed above, LISA functions as a pair of two-arm, Michelson detectors, which we label I and II. The output 
from these two detectors can be represented by the vector Sa{t) (with a = /, //). In what follows it will be convenient 
to work with the Fourier transform of the signal; the convention we use is 

/>oo 
e^^'f's^{t)dt. (34) 
-oo 

Now, the output Sa{t) is the sum of incident gravitational waves ha{t) and instrumental noise na{t). For simplicity 
we assume that (i) the noise is stationary and Gaussian, (ii) the noise in detectors I and II is uncorrelated, and (iii) 
the noise spectral density Sh{f) is the same in the two detectors. 'Gaussianity' means that each Fourier component 
ha{.f) has a Gaussian probability distribution. The combination of assumptions (ii) and (iii) is expressed by the 
relation 

{nM) MfT) = ls{f - f')Sr^{f)d^0, (35) 

where '( )' denotes the "expectation value" and <?„(/) is the (single-sided) noise spectral density for each detector. 
Under the assumptions of stationarity and Gaussianity, we have a natural definition of the inner product on the 

vector space of signals [34]: 

/•oo 

(p|q)^2^ / \p*M)QM)+PM)rM)]/Sn{f)df, (36) 

where Pa{t) and qa{t) are any two signals. One can show, based on this definition, that the inner product of pure 
noise n with any signal is a random variable with zero mean and unit variance. In particular, the probability for the 
noise to have some realization no is just 

p(n = no) oc e-("°l"°)/2^ (37) 

Thus, if the actual incident waveform is h, the probability of measuring a signal s in the detector output is proportional 
to e~(^~''l^~''^/^. Correspondingly, given a measured signal s, the gravitational waveform h that "best fits" the data 
is the one that minimizes the quantity (s — h | s — h) . 

The SNR for an incident waveform h filtered by a perfectly matched template T = h is 

SNR[h] ^^^^] (38) 
rms (hjn) 

where rms (h|n) is the rms value for an ensemble of realizations of the detector noise, n. From the definition (36) it 
follows [34] that for any two signals Pa{t) and qa{t), the expectation value of (p|n)(q|n) is just (pjq). In particular, 
we have rms (h]n) = (h]h)^/^, and hence the SNR of the detection is approximately given by just 

SNR[h] = (h|h)i/^ (39) 

For a given incident gravitational wave, different realizations of the noise will give rise to somewhat different best-fit 
parameters. However, for large SNR, the best-fit parameters will have a Gaussian distribution centered on the correct 
values. Specifically, let A" be the "true" values of the physical parameters, and let A" + SX"" be the best fit parameters 
in the presence of some realization of the noise. Then for large SNR, the parameter-estimation errors (5A° have the 
Gaussian probability distribution 

p{SX") = Afe-i^-"^^''^^'. (40) 
Here Tab is the so-called Fisher information matrix, defined by 



dh 



dX^J' 



(41) 
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and N = ^/det{T/2■K) is the appropriate normalization factor. For large SNR, the variance-covariance matrix is given 

by 

{SX"-SX^) = (r-i)"'' + ^^(SNR)-!. (42) 

1/2 

We define AA"^ = {SX°-SX°-) ' (the repeated index is not being summed here). The uncertainty in the source's angular 
position, AQs (a solid angle), is given by [10]^ 

AQs = 27r ^J{\iisAcl)sY -{SiisHsf ■ (43) 

The second term in brackets in Eq. (43) accounts for the fact that errors in jjLg and (jyg will in general be correlated, 
so that the "error box" on the sky is elliptical in general, not circular. The 2tt factor on the right-hand side (RHS) of 
Eq. (43) is our convention; with this definition, the probability that the source lies outside an (appropriately shaped) 
error ellipse enclosing solid angle Af2 is e~^^/^^s ^j^g game way, the error ellipse for the spin direction, Afi^f, is 
given by 

Afix = 27r ^{A^lKA4>Ky-{5^lKHK?■ (44) 

The actual inner product, Eq. (36), is formulated in the frequency domain. For a white noise [i.e., <?„(/) = constant], 
the inner product is equivalent to 25"^ Pa{t)qa{t)dt, by Parseval's theorem. Motivated by this formula, we 

shall adopt the following approximate version of the inner product in calculating the Fisher matrix: First, we define 

the "noise-weighted" waveform 

h^{t) = ^h^,n{t)/Sl/^Ut)), (45) 

n 

where we take 

U{t)^ni^{t)+^{t)/n. (46) 
Then we approximate the covariance matrix, Eq. (41), as 

Tab = 2 V r dMt)dbK{t)dt . (47) 

a Jo 

That is, we simply re- weight each harmonic by the square root of the inverse spectral density of the noise, and 
thereafter treat the noise as if it were white. 

The decision to set fn{t) = n[u{t) + (2/n)(7(f)/27r)] is something of a compromise: The radial orbital frequency is v, 
the azimuthal orbital frequency is iy{t) + {"/ /2t:), and the fully correct waveforms will contain harmonics of both. Our 
"compromise" approximates the signal as harmonics of (2/n)(7/27r), which lies between the radial and azimuthal 
frequencies^ and is the "correct" choice for circular motion (in which case the n = 2 harmonic dominates the GW 
output, and only the azimuthal piece enters the waveform). 

V. SNR ESTIMATES 

Our analysis of the SNR build-up follows, basically, that of Finn and Thorne [16]. The main advance here is, of 
course, the fact that we consider realistically eccentric orbits, whereas Finn and Thorne confined their analysis to 
circular orbits. Unlike in the circular-orbit case, where the contribution from the n = 2 harmonic always dominates 
the SNR, eccentric orbits have their emitted power (and contribution to SNR) distributed among higher n-harmonics. 
One of the goals of this section is to explore this mode distribution, for realistic values of the orbital eccentricity. 



^Note Eq. (3.7) of [10] is erroneous. However, the results quoted in that paper are based on the correct expression, Eq. (43) 
here. 

^except for the n = 1 harmonic, which, however, contributes very little to the GW signal and to the overall SNR — cf. the 
plots in Figs. (5)-(9) below. 
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A. LISA noise model 



LISA'S noise has three components: instrumental noise, confusion noise from short-period galactic binaries, and 
confusion noise from cxtragalactic binaries. Our treatment of these three sources follows Hughes [35], though we use 
somewhat different estimates for the sizes of the three pieces. 

For LISA'S instrumental noise, we use the following analytic fit by Finn and Thorne [16]'', based on the 

noise budget given in the LISA Pre-Phase A Report [1] : 

= 9.18 X 10-^V"^ + 1-59 X 10"^^ + 9.18 x 10"^^^ Hz-\ (48) 

where the frequency / is to be given in Hz. 

Next we turn to WD confusion noise. Any isotropic background of indistinguishable GW sources represents (for 
the purpose of analyzing other sources) a noise source with spectral density [36] 

Sr'if) = ^r'pcnawif), (49) 

where pc = 3-ffo/(87r) is the critical energy density needed to close the universe (assuming it is matter-dominated) 
and ^GW = {Pc)^^dpGw/d{ln f ) is the energy density in gravitational waves (expressed as a fraction of the closure 
density) per logarithmic frequency interval.^ For the cxtragalactic WD background. Fanner and Phinney [37] estimate 
that, for / near ImHz, nawif) = 3.6 x 10-^^{f /IQ-^Bzp^ [at Hq = 70km/(sec - Mpc)], so 



Sr = 4.2 X 10-4^ (^^) 



-7/3 

Hz-^ (50) 



Note Eq. (50) is not a good fit to S*^ ' ^ for / ^ 10~ Hz, where mergers cause the spectrum to decrease more sharply. 
However, at such high frequencies, instrumental noise dominates the total noise in any case, so for our purposes the 
extrapolation of Eq. (50) to high frequencies is harmless. 

A recent calculation of the galactic confusion background by Nelemans et al. [38] yields an ^^q^ that is 5.0 x 10"'^ 
times larger than fi^jy^"' (near 1 mHz) [37]; therefore^ 

Sf{f) = 2.1 X lO-'*^ j Hz-i. (51) 

This is larger than instrumental noise in the range ~ 10~'*-10~^ Hz. However, at frequencies / ^ 3 x 10"'^ Hz, 
galactic sources are sufficiently sparse, in frequency space, that one expects to be able to "fit them out" of the data. 
An estimate of the resulting noise is [35] 

Sr'+'^\f) = min{sr\f)/eM-'^T-^siondN/df), SrV) + Sf{f)} • (52) 

Here dN/df is the number density of galactic white dwarf binaries per unit GW frequency, Tmission is the LISA mission 

lifetime (so A/ = 1/Tmission is the bin size of the discretely Fourier transformed data), and k is the average number 
of frequency bins that are "lost" (for the purpose of analyzing other sources) when each galactic binary is fitted out 
(k is larger than one because EISA's motion effectively smears the signal from each binary over several frequency 
bins). The factor GWi—nT~l^^i^^dN/df) is therefore the fraction of "uncorrupted" bins, where instrumental noise still 
dominates. For dN/df we adopt the estimate [35] 

dN „ _ _i /iHz^ ^^/^ 



2xlO-=^Hz-i( — ) , (53) 



'^To obtain this expression, we used the expression given at the beginning of page 8 of Ref. [16], where for [ftf ;w6,i!/r(/)]^ 'we 
used the expression given in footnote [44] tiierein, and for A/ we have taken 1/ljyr. 

*Note the RHS of our Eq. (49) is a factor | as large as the RHS in Eq. (3.4) in [36]; this difference arises simply because the 
angle between any two LISA arms is 7r/3 (instead of the n/2 for LIGO's arms), and sin^(7r/3) = 3/4. 

®Note our prefactor 2.f x 10"'*^ is a factor ~ 25 lower than the prefactor cited in Hughes [35], based on his private commu- 
nication with S. Phinney. This large discrepancy seems to be the product of the following two factors. First, it was based on 
the estimate of by Webbink and Han [39], which is 3 times larger than the result of Nelemans et al. [38]. Second, it 

contained a factor 20/3 error due to a misunderstanding of Phinney's normalization convention. 
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and take k.T^^^^^^^^^ = 1.5/yr (corresponding to Tmission ~ Syr and n w 4.5 [40]). To obtain the total LISA noise, we 
just add to Eq. (52) the contribution from the extragalatic confusion background, Eq. (50): 



Shif) = Sr'+^^\f) + '^\f) . (54) 



B. SNR estimates for inspiral orbits 

For a Keplerian orbit, the source is strictly periodic and hence the GWs are at harmonics /„ = nv of the orbital 
frequency. However, as discussed in Sec. IV, to partially compensate for the fact that the radial and azimuthal 
periods are different, we set /n(i) equal to nv{t) + j{t)/7r. Assuming the contributions from the various harmonics 
to be approximately orthogonal, we may approximate the SNR from a single synthetic 2-arm Michelson detector 
(denoted here by SNRi) as 

'^"''■''^' = "-/ (20ffl}!t'.f..) ''''°^'- "'^ 

where "(. . .)sa" means "sky average", i.e., average over all source directions (the factor 20/3 in the denominator 
results from this averaging). The characteristic amplitude hc,n is given, following Finn and Thorne^", by 

hc,n = inD)-'^2EJfn, (56) 
where En is the power radiated to infinity by GWs at frequency /„. To lowest order, this is 



Er. = ^l,'M''/\2n^r/'g{n,e), (57) 





where g{n, e) is given by [6] 



yj4 2 2 

g{n,e) = —{[Jn-2ine) - 2e J„_i(ne) + - J„(ne) + 2e J„+i(ne) - J„+2(ne)] 

+(1 - e2)[J„_2(ne) - 2 J„(ne) + J„+2(ne)]2 + A^[J^{ne)]^} . (58) 

In this section, our major motivation is to investigate the effect of non-zero eccentricity. For this reason, we ignore 
the effect of the MBH's spin — effectively assuming the MBH is Schwarzschild. 

The curves in Figs. 4-9 show the buildup of SNR with time, for each harmonic. As is customary in the LISA 
literature, our plots actually give SNRi the SNR from a single 2-arm Michelson; the actual LISA SNR buildup will 
be a factor ~ times larger. The curves are derived as follows. We use our PN Eqs. (30) and (28) (with S set 
to zero) to evolve e{t) and ^{t) forward in time, up to the point when the CO plunges over the top of the effective 
potential barrier. For a point particle in Schwarzschild, the plunge occurs at amin = M{6 + 2e)(l — e^)~^ [7], so we 
set 

^y^ax = (27rM)-i[(l - e2)/(6 + 2e)f' . (59) 

We denote by e/ the "final" value of e, i.e., the value of e when v reaches the plunge frequency fmax- Then, for each 
harmonic n we use our solution {v{t), e{t)} along with Eqs. (56), (57), and (58) to determine /ic,n(/n)- 

The upper "signal" curves in Figs. (5)-(9) show /ic,n(i) for each n; we "cut off" each curve at fn,max = numax + 
Tr~^^{i'max,^f)- Marks along each curve indicate (from right to left) one, two, five, and ten years before the final 



^"in their SNR estimates, Finn and Thorne [16] tend to consider the quantity fec,m (in their notation), the "modified" 
ciiaracteristic amplitude, introduced to account for the reduction in the GW signal near the plunge, where the available 
bandwidth becomes very small. Here we rather consider hc,m itself: This quantity has the convenient characteristic that when 
integrated against the frequency [through Eq. (55)] it yields the SNR (squared). An estimate of the SNR based on the plots 
given in [16], which show the modified amplitudes rather than hem itself, actually takes into account twice the effect of 
the final plunge: The fact that the frequency changes rapidly near the plunge is already accounted for in the definition of /ic,m, 
just above Eq. (2.2) therein. 
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FIG. 4. GW signal from a WMq CO spiralling into a (non-spinning) 10® Mq MBH at D = 1 Gpc: case where eccentricity at 
the last stable orbit (LSO) is eLSO = 0.3. The curve labeled 'Sinst' shows EISA's sky-averaged instrumental noise level, h^^^^if)- 
The dashed line is an estimate of LISA's overall noise level, h„{f), including the effect of stochastic-background "confusion" 
due to WD binaries (both galactic and extra-galactic). The convex curves show the amplitudes hc^„ of the first 10 n-harmonics 
of the GW signal, over the last 10 years of evolution prior to the final plunge. Along each of these curves we marked 3 dots, 
indicating (from left to right) the GW amplitude 5, 2, and 1 years before the plunge. The orbital eccentricity 10, 5, 2, and 1 
years before plunge is 0.77, 0.67, 0.54, and 0.46, respectively. The orbital frequency 10, 5, 2, and 1 years before plunge is 0.23, 
0.41, 0.70, and 0.94 mHz, respectively. The frequency at the LSO is 1.65 mHz. 



plunge. The lower "noise" curve depicts hn{f) = [{20/3)fSh{f)Y^'^, noise amplitude per logarithmic frequency 

interval [the factor of 20/3 comes from sky-averaging; Finn and Thorne [16] define = {20/3)Sh{f)]. For 

comparison, the instrumental contribution, /i™'(/) = [{20/3) fnS]^^*'{fn)]^^^ , is also plotted. With these conventions, 
the contribution to the SNR from each harmonic is 

(SNRi)2= j{K^n/Kfd{\nf), (60) 

so using the curves one can "integrate by eye" to estimate the contribution to the SNR from each harmonic, and to 
detect when (i.e., how long before the final plunge) most of the contribution is accumulated. 
Here are some points to pay attention to when examining the SNR plots in Figs. 4-9: 

• Systems with MBH's mass of ~ 10^ are ideally "located" in the LISA band. Systems with M = 10^ radiate at 
frequencies where WD confusion noise would likely obscure the capture signal. 

• Systems with M « 10*^ at D = 1 Gpc are detectable with SNRi of 5(/i/1Mq), assuming 1 year of signal inte- 
gration. Combining both synthetic Michelsons, and for two years of integration, LISA's SNR is ~ \Q{ijl/1Mq). 

• As expected, the higher the orbital eccentricity, the more SNR is contributed by high rt-harmonics. At clso = 
0.3, the contribution from n = 3,4 is equally important to that of n = 2. At clso = 0-5j the dominant 
contribution comes at n = 4, and one needs to sum the contributions of at least a dozen modes in order to 
properly estimate the overall SNR. 

• A related point: For 10® M© MBHs and ~ lA/© COs, the last year prior to plunge contributes only a small 
fraction of the potential SNR. Even for ^ lOM© COs, the contribution between 10 and 1 years prior to plunge 
can easily exceed the contribution from the last year. 

The last effect may be further visualized by looking at the total amount of energy radiated in GW up to a time 
t prior to plunge, as a function of t. This is demonstrated in Fig. 10 for a mass ratio of (1OM0)/(1O®M0) and for 
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FIG. 5. Same as in Fig. (4), but for inspiral of a \Mq CO into a 10 M© MBH. The orbital eccentricity 10, 5, 2, and 1 years 
before plunge is 0.46, 0.40, 0.35, and 0.32, respectively. The orbital frequency v 10, 5, 2, and 1 years before plunge is 0.94, 
1.16, 1.39, and 1.51 mHz, respectively. The frequency at the LSO is 1.65 mHz. 
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FIG. 6. Same as in Fig. (4), but for inspiral of a IOM0 CO into a lO^M© MBH. The orbital eccentricity 10, 5, 2, and 1 years 
before plunge is 0.324, 0.313, 0.305, and 0.303, respectively. The orbital frequency 10, 5, 2, and 1 years before plunge is 0.151, 
0.158, 0.162, and 0.164 mHz, respectively. The frequency at the LSO is 0.165 mHz. 
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FIG. 7. GW signal from a lOM© CO spiralling into a (non-spinning) W^Mq MBH: case of a circular orbit (compare with 
Fig. 4 of Finn and Thorne). Notation is the same as in Fig. (4). Note that in our model, a CO in a circular orbit emits GW 
only into the n = 2 harmonic. The orbital frequency 10, 5, 2, and 1 years before plunge is 0.58, 0.74, 0.99, and 1.22 mHz, 
respectively. The frequency at the LSO is 2.20 mHz. 
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FIG. 8. Same as in Fig. (4), except that the LSO eccentricity is taken to be 0.1. The orbital eccentricity 10, 5, 2, and 1 years 
before plunge is 0.38, 0.30, 0.22, and 0.18, respectively. The orbital frequency 10, 5, 2, and 1 years before plunge is 0.49, 0.67, 
0.94, and 1.17 mHz, respectively. The frequency at the LSO is 2.06 mHz. 
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FIG. 9. GW signal from a IOMq CO spiralling into a 10® Mq MBH: case where eccentricity at the last stable orbit is 
clso = 0.5. Notation is the same as in Fig. (4). The orbital eccentricity 10, 5, 2, and 1 years before plunge is 0.88, 0.81, 
0.71, and 0.64, respectively. The orbital frequency 10, 5, 2, and 1 years before plunge is 0.151, 0.265, 0.490, and 0.678 mHz, 
respectively. The frequency at the LSO is 1.13 mHz. 

a variety of LSO eccentricities. For this plot, we used the leading-order expression given in [6] for the total power 
radiated from all n-harmonics: 

E(t) = — (^/M)2[27rz.(i)Af]i"/3[l - e{tf\-'l'^ [l + (73/24)e2(i) + (37/96)e^(i)] . (61) 
o 

The percentage of energy radiated up to time t, out of the total energy radiated during the capture, is then calculated 
through 



%E,Ut) = I 1 - ) X 100. (62) 



The results are striking: For a 1O^M0 MBH and a IOMq CO, about half the total GW energy is released earlier than 
10 years before the final plunge; for a IMq CO, half the energy is emitted already 100 years before plunge. Most 
of this energy is released in the LISA band, in short "spurts" as the CO passes close to the MBH. This provides a 
vivid demonstration of the potential "threat" imposed by self-confusion. A systematic analysis of this issue will be 
presented elsewhere. 

Finally, recall that all of our results regarding SNR values should be regarded as ballpark estimates rather than 
definite, due to the many approximations made here. Here is a summary of the various approximations involved: (i) 
We incorporate an approximate model of LISA's instrumental noise, which is inaccurate at the high-frequency end; 
(ii) We model the source using PN evolution and quadrupole emission — both approximations being worse, again, at 
high frequencies; (iii) We assume a non-rotating MBH; (iv) We assume the contributions from the various n-modes 
to the SNR are uncoupled; and (v) We assume that an ideal, coherent matched filtering search is carried out over the 
entire observation period. Nevertheless, we believe our results amply illustrate the importance of higher harmonics 
and of searches with multi-year integration times. 

C. SNR for a low-mass main-sequence star at Sgr A* 

Freitag [18,41] recently pointed out that LISA might be able to detect a few ~ O.IM0 main-sequence stars captured 
by the MBH at the center of our own galaxy. The strongest such source will still have ^ 10^ years to go before plunge. 
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FIG. 10. Percentage of energy radiated (out of the total energy radiated in GW to infinity during the entire capture process) 
as a function of the time left to the final plunge. The CO and MBH masses are WMq and W^Mq, respectively. This plot was 
generated using Eq. (62) with E approximated using Eq. (61). Note how large the time scale is for energy to be emitted in 
GW by the CO-MBH system, compared to LISA mission time (~ 2 years). 

and so will currently be at the low end of the LISA frequency range, but because it is so close to us, Freitag estimated 
that it could still yield a SNRi as high as ~ 100. We have re-examined this estimate, using the example from Fig. 1 in 
Freitag [41], in which /i = O.OGMq, M = 2.6- IO^Mq, D = 8 kpc, and e = 0.8 at ten million years prior to plunge. The 
results are shown in Fig. 11. Since for the galactic source considered here, there is practically no frequency evolution 
during the observation time, our convention in Fig. 11 differs somewhat from those in Figs. 4-9. As in Figs. 4-9, 
the different curves correspond to the contributions from different harmonics, but here we plot hc^n times the factor 
[{f~^dfn/dt) X 2yr], so the height of the marked point above the noise curve now gives SNR accumulated over 2 years 
of observation. The marked points correspond to 2-yr observations carried out 10^, 5 • 10®, 2 • 10®, 10®, and 10^ years 
before plunge, respectively. (Actually, there is no "plunge," as the low-mass star is tidally disrupted at least ^ 50 yrs 
before plunge would occur, but we can ignore that for this analysis.) 

For a signal coming 10® years before plunge, we find a 2-yr SNRi of ~ 11 (the square root of the sum of the squares 
of the contributions from all modes). A two-year observation only 10^ years before plunge would yield SNRi « 55. 
These SNRi values are a factor ^ 4 times smaller than the values obtainable from Fig. 1 of Frietag [41] (when the 1-yr 
results he gives there are scaled up to two years), but still support the idea that this is a potentially observable source. 
As Freitag points out [41], tidal heating could possibly (under pessimistic assumptions) disrupt the low-mass star 
sometime between 10® and 10^ years prior to plunge, so the higher SNRi value must be treated with some caution. 

VI. PARAMETER ESTIMATION 
A. Numerical implementation 

For our parameter-accuracy estimates we wrote a simple numerical code, based on the following prescription. 

• Pick a specific point A" in parameter space. We found it convenient to first prescribe a value for the eccentricity 
at the LSO, clso; then get the associated LSO frequency using, for simplicity, the Schwarzschild value [7] 

ui^so - (27rM)-i[(l - eiso)/(6 + 2eLso)f^; (63) 

and finally obtain cq and vq by integrating the evolution equations of Sec. Ill F one year back in time. The 
parameter to is set to be in the middle of the integration time. Namely, the integration is carried out from 
Unit = to- (l/2)yr to tLSO = + (l/2)yr. 
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• For that point, calculate the n-modes hj,n{t) and hii,n{t) of Eq. (14), for the last year of inspiral, as described 
in Sec. III. The number of modes one has to take into account varies with the prescribed LSO eccentricity. 
We determined this number by requiring that the relative error in the final Fisher matrix components due to 
omission of higher-n modes is not greater than ~ 10~^. At glso — 0.5 this meant summing over ~ 20 modes. 
Note that at the n'th mode, the 1-year long function h^^n contains a huge number of wave cycles — roughly 
~ lO^n. The time resolution has been set such that each wave-cycle is sampled at least 10 times. At rt = 20, 
this meant a time resolution of about 5 seconds. 

• For each relevant rt, use Eq. (54) to calculate LISA'a noise Sn{fn)- Then sum over modes using Eq. (45) to 
obtain the "noise- weighted" waveforms hi{t) and hij{t). 



Calculate the SNR through Eq. (39), using the time-domain approximation for the inner product: 

■tLSO 



/■ILSO 

SNR2 = 2 ^ / hl{t)dt. 



(64) 



• Calculate the 13 x 2 derivatives daha = dha/dX"". We take these derivatives numerically, through 

KiX" -^X" + 5X''/2) - KiX" ^ X" - 5X''/2) /SX" (for each A'^ = A°, . . . , A^^). Namely, for each derivative 

we calculate the waveform twice, with the relevant parameter shifted by ±^A'^/2. The shift 5A° is set such that 
the resulting relative error in the Fisher-matrix elements is <^ 10~^. [The dependence of the waveform in some 
of the parameters is such the derivatives can, in principle, be taken analytically. For some other parameters 
(like eo or <I>o) the dependence is less explicit, and one is required to take the derivatives numerically. We found 
it convenient (and accurate enough) to take all derivatives numerically.] 

• Calculate all elements of the Fisher matrix Tab, using Eq. (47). 

• Use Eq. (42) to calculate the measurement error for each of the parameters: 

AA" === V(r-i)"" (65) 

(no summation over a on the RHS). Obtain Ails and Afljc using Eqs. (43) and (44). To invert the Fisher 
matrix, we used a numerical subroutine based on dgaussj ( ) of Ref. [42] . 
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• Finally, multiply each of the various AA° by SNR/30. This, in effect, normalizes the distance to the source such 
that the SNR becomes 30 (roughly the SNR output from a 3Mq CO captured by a lO'^M© at £> = 1 Gpc). 

B. Results for 1-yr integrations 

Although in Sec. V we strcssc^d the importance of the signal-to-noisc built up in the last several years prior to 
plunge, in this section, due to limitations of computer memory and speed, we restrict ourselves to waveforms coming 
from the last year of inspiral. In essence, in this section we pretend that LISA is "off-line" prior to tLSO — 1 yr, where 
tLSO is the instant of plunge. 

We presents our results from inverting the Fisher matrix for MBH mass of lO^M©, CO masses of IM© and lOM©, 
and a range of values of the MBH spin and the orbital eccentricity at the LSO. The rest of the parameters have been 
set as follows: to = Ilso — (l/2)yr (middle of integration), 70 = 0, $0 = 0, 9s = 7r/4, 05 = 0, A = 7r/6, ccq = 0, 
Ok = 7r/8, and (1)k = 0. Again, the angles (po and cto specifying LISA's position and orientation at to are set to 
zero. Tables II and III give the results. We expect that the measiircmcnt accuracies for intrinsic parameters will not 
depend very sensitively on the source's position and orientation, but that the measurement accuracies for extrinsic 
parameters (e.g., LISA's angular resolution) will depend rather more sensitively on the actual values of these angles. 
(This was found to be the case in LISA measurements of MBH-MBH coalescences [21], and we have also been able to 
verify it, to a limited extent, in the LISA capture case.) 
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TABLE II. Parameter accuracy estimates for inspiral of a IMq CO onto a 10® Mq MBH at SNR=30 (based on data collected during 

the last year of inspiral). Shown are results for various values of the MBH's spin magnitude 3 and the final eccentricity blso- The rest 
of the parameters are set as follows: to = tLSO ~ (l/2)yr (middle of integration), 70 = 0, <J>o = 0, 6s = 7r/4, tpg = 0, X = 7r/6, ao = 0, 
Ok = n/8, 4>k = 0. 

Summarizing the results in Table III, for the case of a lOM© BH spiraling into a IO^Mq MBH, the measurement 
accuracies for the seven intrinsic parameters are (very roughly) A(lnM) ^ 2xl0~^, A(ln/i) ~ 10~*, A(S'/M^) ^ 10~^, 
A(cosA) - 10-3, Aeo ~ IQ-^, and A$o ~ A70 ~ 0.5. 

Comparing Tables II and III {IMq vs. lOM© CO), we see A(lnM), A(ln/i), A{S/AP), A(cosA), and Aeo are all 
typically about an order of magnitude smaller for the IOMq CO case (again, for fixed SNR and fixed 1-yr integration 
time). This general trend of better accuracy for higher CO mass is not hard to understand. First, since i> oc /i, it is 
clear that Fin^ in^ should scale roughly like /x^, so A(ln/Lt) = [(F-^)^"** in;iji/2 g]jould scale roughly as /i"^. Similarly, 
the derivative 9inM^a(i) has size of order ha(t)din M^{t) (during most of the integration time), and it is clear from 
Eqs. (27)-(28) that the typical size of dinM^it) scales roughly like /j.. The magnitudes of ds/M^ha{t), dcosxha{t), and 
de„ha{t) also scale roughly linearly with /x, for the same reason. Hence, it is reasonable to expect that errors in these 
five intrinsic variables should scale roughly like fi~^. 

We next turn to the extrinsic parameters. Our few examples suggest that LISA's angular resolution for capture 
sources is Aflg ^ 10^^ radians, while the MBH spin direction can be determined to within AUk ~ 5 x 10"^. We 
find A[ln{n/D)] ~ 5 x lO'^, typically. Since A(ln/y,) < A[ln(/x/Z))], it is clear that A(ln£)) ~ 5 x lO'^ too. [We may 
verify the last statement by examining the [In /j,, ln(/j,/D)] minor of F~^: We have 

/-p—l\\nD,\nD /-p—l\\n{fi/ D) ,\n{fi/ D) , f-p—l\\nfi,\nfi ^/p — l\ln/j,ln(/j/D) 
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TABLE III. Parameter extraction accuracy for inspiral of a IOMq CO onto a 10® Mq MBH at SNR=30 (based on data collected 

during the last year of inspiral). Shown are results for various values of the MBH's spin magnitude S and the final eccentricity clso- The 
rest of the parameters are set as follows: to = tj^so ~ (l/2)yr (middle of integration), 70 = 0, #0 = 0, ds = ''r/4, = 0, A = 7r/6, ao = 0, 
Ok = n/8, 4>k = 0. 

and since typical values are (r-iyn{t,/D)M(M/D) _ ;^q-4^ (p-i^nA^Jn^ _ iq-^^ .^^^ (p-iynAi,in(M/i5) _ 

that indeed A(ln£)) « A[\n{ij/D)].] Finally, as a check, it is easy to see that A[ln(/x/JD)] must be greater than 

SNR~^ = 0.033 (since the signal amplitude is linear in i^/D), which is indeed satisfied in every column of our tables. 

C. Comparison with other results in the literature 

Our angular resolution results can be compared to results by Cutler and Vecchio [43] on LISA's angular resolution 
for monochromatic sources. For a monochromatic source with fgyj = 3 mHz and SNR=30, LISA's angular resolution 
is typically Afi,, ^ 5 x 10^'' (estimated by interpolating between Figs. 2 and 3 in [43], after rcscaling those figures to 
SNR=30), which is only a factor ~ 2 smaller than our result for capture sources. Since capture sources have twice as 
many unknown parameters as monochromatic sources (14 versus 7), it is clear that LISA's angular resolution must 
be worse for the former (at the same SNR), but the "good news" is that this degradation appears to be quite modest, 
based on our limited sample. 

Our results on the mass and spin determination accuracy can be compared to previous results by Ryan [2] and 
Poisson [22]. Ryan's waveforms are based on PN evolution equations (similar to ours), while Poisson's are based on 
a Taylor expansion of the waveform phase near plunge, with expansion coefficients obtained from numerical solution 
of the Tcukolsky equation. Both these authors consider only circular, equatorial orbits (so Cq = and cos A = 1.0, 
a priori). Both simplify the calculation further by ignoring the waveform modulation caused by LISA's motion (so 
they effectively pretend LISA is fixed at the center of our solar system), and by restricting attention to the waveform 
generated by just a single pair of LISA's arms. Thus, their waveforms are determined by only 5 parameters: an overall 
amplitude and phase, the two masses, and S/M^. (Clearly, these simplifications were intended to make the Fisher 
matrix calculation essentially identical to the corresponding calculation for LIGO measurements of binary black hole 
coalescences.) 

The fact that Ryan and Poisson effectively "toss out" most of the unknown parameters obviously tends to decrease 
the calculated error bars for the included parameters. On the other hand, their highly simplified waveforms obviously 
carry much less information than the true waveforms, which tends to have the opposite effect. A priori, it would seem 
difficult to guess whether the net effect of their approximations is to underestimate or overestimate A(lnM), A(ln^), 
and A{S/M'^). Therefore, unfortunately, their work does not seem to provide a useful check on ours. Nevertheless, 
Ryan's and Poisson's papers were an interesting first-cut at the parameter estimation problem, and it seems interesting 
to compare our results to theirs. 

For a IOM0 CO and IO^Mq MBH, Ryan [2] obtains (at SNR=30): A(lnM) = 1.8 x 10"^, A(ln/i) = 1.9 x 10-^ 
and A{S/M'^) = 4.9 x 10"''. For the same masses and SNR, Poisson [22] states the resuhs A(lnM) = 6.7 x 10"^, 
A(5/M^) = 1.7 X 10~^, and A(ln?7) = 1.8 x 10"'^, where ry = ^/M. Since both Ryan and Poisson consider only the 
case eo = and evaluate the Fisher matrix at the point S/M"^ = 0, their results are most usefully compared to those 
in column 1 of our Tables (i.e., eo = 0.1 and S/M^ = 0.1). Our estimates of mass-determination accuracy are within 
roughly an order of magnitude of those quoted by Ryan and Poisson, and in fact lie between them. [Poisson does not 
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quote a result for A(ln;u), but we can still compare our results directly to his by using 



^p— iyn?7,ln77 _ ^p — 1 /t,ln /:i _j_ ^p— iynM,lnM 2^p~l /^il" ^ (66) 

which gives Ar] — 3.0 x 10~^ for our case] Our A(5/A'/^) is ^ 10 times smaller than Ryan's and ^ 50 times smaller 
than Poisson's. We guess this is because in Ryan's and Poisson's waveforms the MBH spin affects only the orbital 
phase (whereas in the true waveforms the spin also controls the Lense-Thirring precession rate), leading them to 
overestimate the covariance of the spin with the two mass parameters, and hence to overestimate A{S/M^). 



D. Parameter extraction for a low-mass main-sequence star at Sgr A*. 

Finally, it's interesting to consider the parameter extraction accuracy for captures of LMMSs at the center of the 
Milky Way — the type of source whose anticipated SNR we discussed at the end of Sec. V. In Table IV we consider 
a few possible low-frequency orbits of a O.O6M0 CO around the 2.6 • W^Mq MBH at the known distance and sky 
location of Sgr A*. We take D=7.9 kpc [44] and 9s = 1.66749; the value for (ps is picked arbitrarily, since the orbital 
location of the detector during the observation is unknown. The columns with eccentricity 0.43 and frequency 0.0035 
mHz correspond to the orbit whose SNR output has been discussed above (in Fig. 11; also [41]), and assume the 
observation is made 10^ years prior to plunge (just before tidal forces play an important role). In other columns we 
consider other possible values of the eccentricity, frequency, and MBH's spin. In calculating the Fisher matrix, we 
assumed a data integration time of 2 years. 
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TABLE IV. Parameter extraction accuracy for a low-mass main-sequence star at Sgr A*. We assume M = 2.6 • 10® Mq, // = 0.06Mq, 
and data integration lasting 2 years. We also assume the star is observed a million years before the (theoretical) plunge, just before tidal 
effects bccoiuc important. Each column of the table refers to a different choice of the MBH's spin, orbital eccentricity e and frequency 
u at the time of observation. The other parameters are set as follows: 70 = 0, 3>o = 0, &s = 1-66749 (true value for Sgr A*), <ps = 0, 
A = 7r/6, ao = 0, Ok = '""/S, ipK = 0- Most of the values given in the table result from inverting the full, 14 x 14-d Fisher matrix. The 
values for AQk obtained by inverting the 11 X 11 minor that excludes the CO's mass fi and the two sky-location coordinates 9$ and <ps 
(whose precise values are known for Sgr A*) are given in parentheses. (For all other parameters, using the known sky position did not 
signficantly improve measurement accuracy.) 

Most of the values in the table result from inverting the full, 14-d Fisher matrix. We find that the sky location 
of the source can be determined to within ~ 0.01 steradians. Unfortunately, the distance D to the source is entirely 
degenerate with the CO's mass, and the latter cannot be determined by itself from the GW signal. [This is as expected, 
since the accumulated effect of radiation reaction on the waveform phase, after observation time Tobs ~ 1 yr, is merely 
A(j) ^ ''^fT^bs ^ '^■f'^obs/ i^^^y^) ^ 10^^, for / ^ 10"'* Hz.] Hence, the distance to the galactic center cannot be 
determined by the GW signal alone. However, in case where a detected GW signal appears to come from a capture 
at Sgr A* (as its sky position is consistent with the galactic center and its / is extremely small), one could use the 
value of D known from astronomical measurements, in order to determine the CO's mass /i. (Recent studies of stellar 
dynamics and RR Lyrae stars have specified D to within 4% [44]). The relative error in /j. will then be approximately 
A[ln(/i/£))], or roughly 10%. 
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Once a source is confirmed to be at the galactic center, we may eliminate the sky location Qs (known exactly for Sgr 
A*) from the search, in order to improve the accuracy in determining the other parameters. We may also eliminate 
the CO's mass /U from the parameter list (while keeping i-i/D) — again, because radiation reaction is negligible for this 
source. 

To estimate the effect of this extra information on parameter extraction, we inverted the 11-d minor of the Fisher 
matrix obtained by excluding the rows and columns associated with /zg, (j)s, and In/i. We found that this produces 
only a negligible improvement in measurement accuracy for the rest of the parameters, except for the MBH spin's 
direction; AQk decreases by a factor ~ 1.5-3. These improved values of AQk are given in the table as well, set off 

by parentheses. 

The more interesting results of Table IV concern the intrinsic parameters of the galactic hole. We estimate that the 
mass of the central MBH (as well as the orbital eccentricity) could be determined to within relative error ~ 4 x 10~^ 
to 2 X 10""'^, i.e., comparable to or better than what is achieved today using astromctric methods. More impressive is 
the determination of the MBH's spin from the GW waveform to within ~ 0.005 (in units of M^). It is hard for us to 
imagine an alternative method that would allow such an accurate measurement of the spin of the MBH at the Milky 
Way's center. 

Finally, we point out that two approximations we have made — using PN equations of motion and a low-frequency 
approximation to LISA's response — should be really quite accurate for LMMS captures in our galaxy, since these 
inspirals are viewed at a very early stage (the orbits considered in Table IV have their pericenters at ^ 20M), and 
most of the SNR comes from frequencies in the range 10~''-10~^ Hz. 

VII. CONCLUSIONS, CAVEATS, AND FUTURE WORK 

Realistic capture orbits will be nonequatorial and somewhat eccentric, in general. We have made a first cut at 
answering some LISA data analysis questions, for such realistic cases. 

The figures in Sec. V illustrate how the LISA SNR builds up over time for eccentric orbits (assuming there, for 
simplicity, that the MBH is nonspinning) . These figures show that, for ~ lOM© COs captured by ~ lO^M© MBHs, 
the last few years of inspiral can all contribute significantly to the SNR; for ~ IMq COs, the last few decades can 
be significant. Clearly, these long waveforms will increase the computational burden on matched- filtering detection 
schemes, and will exacerbate the self-confusion problem, since eccentric-orbit inspirals can deliver signicant GW energy 
into the LISA band many years before they become individually detectable. 

The tables in Sec. VI represent our attempt to estimate LISA's parameter estimation accuracy, for captures. For a 
typical source (a IOMq CO captured by a IO^Mq MBH at SNR of 30) we find that A(lnM), A(ln/i), A(S'/M2), and 
Aeo will all be roughly ^ 10~^, while Ailg ~ 10""^ and AUk ~ 5 x 10~^. Due to computational limitations, those 
results are based on measurements from only the final year prior to plunge. We naturally expect the measurement 
errors to decrease when one considers waveforms lasting the entire length of the LISA mission. In this sense, the 
numbers above represent rough upper limits to LISA's measurement accuracy. On the other hand, these estimates are 
based on a small sample of hand-picked points in parameter space. In the future we will improve these estimates by 
doing several- year integrations and a full Monte Carlo sampling of parameter space. Of course, our PN waveforms are 
probably not very accurate for COs very close to the plunge, but still the above estimates arc the best ones available. 

LISA has a reasonable chance of detecting LMMSs captured by the MBH at the center of our own galaxy. In this 
case, the sources will be detected ~ 10^ yrs before the final plunge, when our PN waveforms and our low-frequency 
approximation to LISA's response function should both be more reliable. We find that, in a two-year integration, LISA 
could determine the magnitude and direction of our MBH's spin to within A{S/M'^) ~ 5 x 10"'' and AnK ~ 10"^ 
and measure the mass of the infalling star to within 10%. 
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APPENDIX A: LENSE-THIRRING COMPONENT OF THE PERICENTER PRECESSION 

At IPN order, the piece of 7 = d^y/dt proportional to S is [30] 

= (L-ricosA-5-n)d-87ri^(5/M2)cosA(27rMzy)(l-e2)-3/2, (Al) 

I — {L ■ n)2 

where d = da/dt is given in Eq. (31). [To derive this equataion, use Eq. (4.4.45) of Brumberg [30], together with his 
Eqs. (4.4.41) and (4.4.33), and along with the definitions in Eqs. (1.1.6)-(1. 10.10) therein.] 

The expression (Al) is somewhat complicated (and worse, badly behaved as L approaches n) due to the usual 
convention of defining 7 as the angle from x = [L{L ■ n) — n]/{l — {L ■ fi)'^)^^'^ to pericenter. In the terminology of 
BCV [24], 7 is neither wholly an intrinsic variable, nor an extrinsic one; it is of a mixed type. As explained in Sec. 
Ill, it is preferrable to introduce a pericenter angle that is defined purely intrinsically (i.e., without reference to the 
observer). Accordingly, we define 7 as the angle from L x S to the direction of pericenter. Then 7 and 7 are related 
by 

7 = 7 + /3, (A2) 

where /3 is given by Eq. (21). Not only is 7 wholly intrinsic, but we claim it also obeys a simpler (and better behaved) 
evolution equation, Eq. (29): 

^\^s = -12-kucos\{S/M'^){2'kMu){1 - e^)-^/^ (A3) 

The purpose of this Appendix is to prove this claim; i.e., we show that Brumberg's Eq. (Al) is indeed equivalent 

to our Eqs. (A3) and (A2). 

Since it is obvious from the definitions that 

70 = 70 + /3o , (A4) 

what we need to prove is that 

/3 = 7-4- (A5) 
Prom Eqs. (Al), (A3), and (31), and noting that 7 and 7 differ only in their piece cx S, the RHS of Eq. (A5) becomes 

7 — 7 = {L ■ h cos A — 5 • n) d + cos A d 

1 — (L • n)2 

^ cos X - {L-h){S ■ h) ^ 
1 - (i • n)2 

The left-hand side (LHS) of Eq. (A5) can be written as 

P = {cos(3)-^d{smf3)/dt, (A7) 
which, substituting for sin/3 and cos/3 from Eq. (21), yields 

^ _ [1 - ■ »)-] [ ^■"■^ A [L ■ n) ^ (cos XL ■ // - 5 ■ n)[L ■ n)[L ■ n) \ 
n-{SxL) \ [l-(i.n)2]^/' [l - (L ■ h^f j' 

Using now the evolution equation for L, L = d S* x L, the last expression reduces to the RHS of Eq. (A6). 
This proves Eq. (A5), and hence Eq. (A3) is proved. 
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Wc note Eq. (A3) can also be seen as a special case of Eq. (Al), obtained by setting n equal to S in the latter. 
Finally, we also note that the term 



(1/ • ncos A - 5 • n) d (A9) 

1 — (L • n)2 

in Eq. (Al) for 7 is precisely the time derivative of what has been referred to in Ref. [23] (in the context of quasi- 
circular orbits) as the "Thomas precession phase". [The term (A9) is precisely the RHS of Eq. (28) in [23].] Wc 
see here that this "Thomas precession" term is really just the "L x S piece" of the pericenter precession [while the 
remainder of Eq. (Al) is the "L • S piece"]. 



APPENDIX B: RELATIVE MAGNITUDES OF PN CONTRIBUTIONS 

Our goal here is to gauge the suitability of the PN formulae we use in evolving the CO's orbit within our model [Eqs. 
(27) (30)]. The part of the inspiral relevant for detection by LISA takes place entirely within the highly relativistic 
region right next to the horizon, a regime where the validity of the PN expansion is normally to be suspected. 
Accordingly, we do not expect our PN expressions to show any convergence; we do wish to make sure, though, that 
truncating the PN expansion as in Eqs. (27)-(30) does not lead to a pathological, unphysical evolution within the 
range of parameters relevant for our analysis. 

In Fig. (12) we compare the magnitudes of the contributions from the various PN expansion terms to each of the 
time derivatives 1), e, and 7. We consider the "worse case", in which the CO has just reached the last stable orbit 
(LSO) and is about to plunge into the MBH. For simplicity, we take {2ttMi')lso = [(1 — e^)/(6 + 2e)] , which is the 
value for a non-rotating MBH. 

As expected, the various PN terms are of comparable magnitudes. Nevertheless, our evolution cqiiations do not 
develop pathologies, so long as the LSO eccentricity is not too high. In particular, the orbital frequency increases 
monotonically (no "outspiral"!) and the eccentricity decreases monotonically (the orbit circularizes) throughout the 
entire evolution. 



APPENDIX C: ESTIMATING THE EFFECT OF THE CO'S SPIN 

In this appendix we refer to the spin of the MBH as 5*1, and denote by S2 the spin of the CO (in the rest of 
the paper the MBH's spin is denoted as simply S). We present simple estimates of the effects of S2 on the orbital 
evolution, showing that, at any instant, these are smaller than the effects of Si by an amount of order {h/M){x2/xi)j 
where xi = [^il/M^ and X2 = |<S2|/m^- This implies that over the course of the last ~ lyr of inspiral, the presence 
of the CO's spin modifies the accumulated orbital phase and precession angles by a few radians at most. 

Clearly xi < 1- Theoretical upper limits on X2 are: X2 < 1 when the CO is a BH; X2 <; 2 for a 1.4Mq NS; and 
X2 ^ 5 for a uniformly rotating, 0.5Mq WD [45]. We do not know the spin of any stellar-mass BH, but the NSs and 
WDs we find in nature are typically spinning at well below these maximum possible rates. Most known pulsars have 
X2 ~ 10~^ — 10~^ (recycled pulsars being the major exception), and it appears that WDs typically have X2 ~ 10~^ 
[46]. 

Our estimates here are based on low-order, post-Newtonian equations for circular-orbit binaries. Of course, the 
capture orbits of interest to us can be moderately eccentric in the LISA band, and they are sufficiently close to the 
MBH that higher-order post-Newtonian terms can be comparable in size to the lower-order terms that we do include. 
We expect that these higher-order PN terms, and modifications for e ^ 0, would modify our estimates on the effects 
of S2 by factors of order unity, but we anticipate that our PN estimates at least give the order of magnitude. 

Two other recent works that consider the effect of the CO's spin on the capture waveform are by Hartl [47] and 
Burko [48]. 



1. Effect of L • 5*2 and Si ■ S2 terms on the waveform phase 

For circular orbits, let oj = -I- 7 be the angular frequency of the CO around the MBH. Defining ^ to be the 
accumulated orbital phase from time ti to tf, 

^ = / ujdt= -duj, (CI) 

Jti Juii ^ 
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FIG. 12. Relative magnitudes of the various PN terms. The plots compare the contributions to u (upper left panel), e 
(upper right panel), and 7 (bottom panel) from the various terms in the PN expressions (28)-(30), when the CO is at its last 
stable orbit. The various terms are of comparable magnitude, as expected in this "worse case". Note, though, that so long as 
eLSO < 0.6, the evolution of the orbital parameters at this order is free of potential pathologies [like "outspiral" (;> < 0), or a 
negative pericenter advance (7 < 0)]. 



Kidder [49] finds that the change in the orbital phase 5* due to non-zero 52 is 



1536 



(C2) 



[Actually, Eq. (C2) refers to the case where the dot products L ■ §2, L ■ Si, and 5*1 • 52 are constants (e.g., the case 
where these three angular momenta are all perfectly aligned or anti-aligned, so there is no Lense-Thirring precession). 
For the realistic case, where L and 6*2 undergo significant Lense-Thirring precession, these dot products would be 
replaced by appropriately weighted time-averages.] For our case of inspiral orbits near plunge and (tf — ti) of order 
a year, we have [(Mwj)'^/^ - [Mujfy^/^] - 1 - 2 and [{Muoi)-'^/^ - (Ma'/)^2/3] ^ 3 _ 10. Therefore, when X2 and 
Xi are of order one, A^* is at most a few radians, while 5* itself is of order 10^. 



28 



2. Effect of S2 on pericenter precession, d'y/dt 

The 52 contribution to d'y/dt is given (to lowest nontrivial PN order) by Brumberg [30], in his Eqs. (4.4.41) and 
(4.4.45). We do not reproduce those equations here, but just note that the relevant terms are (for low eccentricity) 
of order X2(M/Af)(Af/r)'^/^a; and xiX2(M/-^)(Ai"/r)^a;, respectively. Therefore, the integrated effect of ^2 on 7 is 
A7 ~ X2if^/M)Norb radians, where Norb is the number of orbits during the integration time. Hence, A7 is of order 
one radian for X2 ~ 1 and typical values n/M 10~^ and Norb ~ 10^. 



3. Effect of S2 on Si and L 



We next estimate how ^2 affects the Lense-Thirring precession. Our estimates are based on the following post- 
Newtonian processional oqiiations given by Apostolatos et al. [23], whic;h wore obtained by averaging the spin-orbit, 
spin-spin, and radiation reaction torques (all calculated to lowest nontrivial post-Newtonian order) over one complete 
orbit: 



ZjJL 



X L - ^1 [(52 • L)5i + (5i • 1)82] X 



5 r \ r / 



^ 2 ^ ^ 1 
Si = —L X 5i + — 



\S2-\{S2-L)L 



X Si 



S2 = ir%—L xS2 + \ 



\si-\{Si-L)L 



X 52 



(C3) 

(C4) 
(C5) 



[Actually, Eqs. (C3)-(C5) are a simplified version of Eqs. (lla)-(llc) in Apostolatos et al, differing by fractional 
corrections of order n/M.] 

Two results that follow from Eqs. (C3)-(C5) are: (i) the two spin magnitudes, |5i| and |52|, are constants of the 
motion; and (ii) the total angular momentum, J = 5i -|- L -|- 52, is also constant, except for the orbital momentum 
that is radiated away: 



5 r \ r J 



(C6) 



Therefore, A J — the total change in J over the observed inspiral — has magnitude of order |AL| ~ {iJ-M). 
Combining this with the fact that the magnitudes of Si , L, and 52 are roughly in the ratio 



|5i| : |L| 



XI : m/M : X2{l^/Mf 



(C7) 



(i.e., |5i| >> |L| >> |52|), one easily sees that the yearly change in the direction of Si can be of order xT^Cm/^) 
radians at most, no matter how 52 and L evolve. Thus our approximation of treating Si as fixed, used throughout, 
remains valid when we include the effects of 52. 

Finally, we consider the motion of L. For 52 = 0, we have seen that L simply processes around Si at rate a given 
by Eq. (31). Because the terms in Eq. (C3) involving 52 are all smaller than the dominant {2/r^)Si x L term by an 
amount of order X2(A*/Af) or (X2/xi)(a*/-^)) they represent a small perturbation on the simple precession picture, 
which can clearly be absorbed into a time-varying A (the angle between Si and L) and a perturbed precession rate 
a. These extra terms in L are of order X2{p/M){M /rY^'^uj and xiX2(M/Af )(Af/r)^ti), respectively — the same as for 
the 52 terms in d'y/dt — and so again integrate up to a yearly difference of order one radian at most. 



4. Conclusion 



If the CO's angular momentum is close to maximal (i.e., within a factor of a few), then the CO's spin is marginally 
relevant for the dynamics over timescales of order a year. However, if the CO is not rapidly rotating, or if one 
is just searching for short stretches (lasting ~ 2 weciks) of the waveforms in the data (e.g., as the first stc;p of a 
hierarchical search), then the CO's spin can be safely neglected. (For such short stretches, in addition to the fact that 
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the waveform phase errors from neglecting ^2 are much less than one radian, it seems likely that these errors can be 
partially compensated for by errors in the other physical parameters.) 
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FIG. 13. Sample equatorial orbits: Shown (left to right, respectively) are 20 minutes, 1 hour, and 6 hours-long samples of 
the go's trajectory just before approaching the LSO. {x,y) is a Gartesian coordinate system in the orbital (equatorial) plane, 
centered at the MBH. The axes give the distance in units of the MBH's mass M. This sequence of plots shows a case with a 
very small LSO eccentricity, e — 0.05. The other physical parameters are set as follows: CO's mass: /i = lOM©; MBH's mass: 
M = 10^ Mq; MBH's spin magnitude: S = M^; 
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FIG. 14. Same as in Fig. (13), for a trajectory with LSO eccentricity e — 0.3. 
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FIG. 16. Sample waveforms. Shown is the responce function hi{t) (as defined in the text) during the last minutes before 
the final plunge. The three panels show cases with LSO eccentricity 0, 0.3, and 0.5 (top to bottom, respectively). The other 
physical parameters are set as follows: CO's mass: fj, = IOM0; MBH's mass: M = IO'^Mq; MBH's spin magnitude: S = M^; 
Angle between MBH's spin and orbital angular momentum: A = 30°. 
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